Flavor ordering of elliptic flows at high transverse momentum 
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Based on the quark coalescence model for the parton-to-hadron phase transition in ultra- 
relativistic heavy ion collisions, we relate the elliptic flow (1)2) of high pt hadrons to that of high px 
quarks. For high pt hadrons produced from an isospin symmetric and quark-antiquark symmetric 
partonic matter, magnitudes of their elliptic flows follow a flavor ordering as («2,tt = V2,n) > (v2,A = 
«2,e) > V2,k > V2,b > («2,0 = v 2,n) if strange quarks have a smaller elliptic flow than light quarks. 
The elliptic flows of high pt hadrons further follow a simple quark counting rule if strange quarks 
and light quarks have same high pt spectrum and coalescence probability. 

PACS numbers: 25.75.Ld, 25.75.-q, 24.10.Lx 



Elliptic flow in heavy ion collisions is a measure of 
the azimuthal asymmetry of particle momentum distri- 
butions in the plane perpendicular to the beam direc- 
tion. It results from the initial spatial asymmetry in the 
transverse plane in non-central collisions and is thus sen- 
sitive to the properties of the dense matter formed dur- 
ing the initial stage of heavy ion collisions . There 
have been extensive experimental P-14J and theoretical 
EKlUJI Hi studies of elliptic flow in heavy ion collisions 
at various energies. For heavy ion collisions at the Rela- 
tivistic Heavy Ion Collider (RHIC), the elliptic flow has 
been measured as functions of the centrality of collisions 
[ 1 1 IH I, as well as the particle transverse momentum 
[Ft ,[^|l4) and pseudo-rapidity jl3). Theoretical studies 
indicate that these experimental results provide not only 
information on the equation of state of nuclear matter at 
high density and temperature |^-^| but also on the scat- 
tering cross section of partons produced in the collisions 



The elliptic flow has also been measured at RHIC for 
different hadron species, such as pions, kaons, nucleons 
and A |2(|^lJ . The experimental data show that at low pt 
the elliptic flow of heavier particles is smaller than that of 
lighter particles. In the hydrodynamical model, this mass 
ordering of elliptic flow at low pt is attributed to the mass 
dependence of radial flow H. For high pt hadrons, we 
expect the flavor dependence to be different from that at 
low pt, since high transverse momentum hadrons origi- 
nate from hard processes while low transverse momentum 
particles are mostly produced from soft non-perturbative 
processes and are much closer to thermal equilibrium. In- 
deed, the observed saturation of hadron elliptic flow at 
Pt > 2 GeV/c pl| , p2| contradicts the predictions from 
the hydrodynamical model §| , but is roughly consistent 
with the results expected from a large parton transport 
opacity jL7| or energy loss (2^] in the partonic matter. 

In this Letter, we shall study the flavor dependence of 
the elliptic flow of high pt hadrons in ultra-relativistic 
heavy ion collisions, using a quark coalescence model to 
describe the phase transition from the partonic matter to 
the hadronic matter. In this model, the elliptic flow of 
high pt hadrons can be expressed in terms of the elliptic 



flow of high pt quarks. As a result, several relations be- 
tween the elliptic flow of hadrons of different flavors are 
obtained. We further discuss some special cases where 
these relations become more transparent. Throughout 
this study, we limit the discussions to hadrons made of 
u, d, s quarks and antiquarks. 

In the quark coalescence model, one assumes that 
quarks and antiquarks are the effective degrees of freedom 
in the parton phase near the phase transition, and they 
combine to form hadrons according to the valence quark 
structure of hadrons. A meson is thus formed from the 
coalescence of a quark and an antiquark, while a baryon 
is due to the coalescence of three quarks. The idea of 
quark coalescence has been used in models such as the 
ALCOR @ or MICOR model || to describe hadron 
abundance and the AMPT model with string melting 
IH to describe the elliptic flow at RHIC. 

In ultra-relativistic heavy ion collisions, high pt par- 
tons are produced from initial hard scatterings between 
nucleons for which the perturbative QCD is applicable. 
From the leading-order calculation, the parton trans- 
verse momentum spectrum from the subprocess of a two- 
parton hard scattering is given by: 



da 1 
— oc 
dt p T 



(1) 



The high pt parton spectrum thus follows an inverse 
power law modulo the corrections from the parton distri- 
bution function in the nucleus and higher-order effects. 
On the other hand, low pt partons, that are produced 
from initial soft processes and dominate the dynamics 
of partonic evolution in heavy ion collisions at RHIC, 
typically have an exponential spectrum close to a ther- 
mal distribution. The parton pt spectrum can thus be 
represented by an exponential function below a certain 
momentum scale po and an inverse power law above pq. 

Let us consider via the quark coalescence model the 
formation of a high pt meson with transverse momen- 
tum ph from one parton with pn and one parton with 
zero pt, or from two partons with equal high pt olpu /2. 
The ratio of the probabilities for forming a high pt me- 
son in these two cases is then proportional to (epr/4po) , 
where n represents the exponent of the inverse power law 
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for final high px partons. Since this ratio is much greater 
than one for px 3> po, a high px meson is dominantly 
formed from the coalescence of one high px parton and 
one soft parton. Similarly, a high px baryon is mainly 
formed from the coalescence of one high px parton and 
two soft partons. 

The transverse momentum distribution F(pr) = 
dN I (dp x dp y ) of initial high px mesons formed after the 
phase transition can thus be expressed in terms of that 
of final high px partons as 



F h { P t) = Fi(p T )cj + Fj(pr)ci 



(2) 



where i and j denote the flavor of the valence quark and 
antiquark of meson H. The coefficient Cj represents the 
capture probability for a soft parton i by a high px parton 
to form a high px meson; it is thus related to the den- 
sity of soft quarks near the phase transition. For high 
px baryons or antibaryons, one can write down a similar 
expression, involving the product of two Cj's, for their 
transverse momentum distribution. 

The elliptic flow is generated during the early stage of 
heavy ion collisions when the pressure gradient and the 
spatial azimuthal asymmetry are the largest P,|5|, |l6| , |l9| . 
In transport model studies, it has been found that the 
elliptic flow in heavy ion collisions at RHIC develops 
mostly in the initial partonic phase, with later hadronic 
interactions having negligible effects on its final value 
HPI . We expect that the elliptic flow of high px hadrons 
are even less affected by hadronic interactions, as the 
proper formation time from a high px parton to a hadron 
is increased by a large Lorentz boost factor in the lab- 
oratory frame, leading to a much lower hadronic den- 
sity when high px hadrons are formed. We can thus use 
Eq. (||) to relate the final elliptic flow of high px hadrons 
to that of high px partons P6|. For mesons, we have 



J cos(20')F H (pi)d(b' 

J F H ( P T)dV 
V2,i(pr)fi(pT)Cj 



V 2 ,j{PT)fj{pT)c l 



Mpt)cj + fj{pi)ci 



(3) 



where 4>' is the azimuthal angle with respect to the re- 
action plane, and /(px) = dN / (2np^dp^) denotes the 
transverse momentum distribution after averaging over 
the azimuthal angle. In the following, we omit the label 
px in the variables «2(px) and /(px) but keep in mind 
that they are evaluated at a given high px . 

For SU(3) hadrons consisting of u, d, s quarks and 
antiquarks, their vi values at high px are then given by: 



«2,7r + : 
V2,4> = 



Vl,ufuC d + V 2d f d C u V2.ufv.Cg + V2,sfsC u 



V2 



P 



fuC d + fd c u 
V 2 ,sfsCs + V 2 ,sf§C s 

fs^s ~l~ fsCs 

fuCd + v 2 ,dfdC u /2 

fuCd + fdCu/2 



fv.Cs ~t~ fsCv 



«2,A = «2,£° 



V2,ufuC d C s + V2,dfdC u C s + V2,s fsC u Cd 
fvC d C s Jd^uCg -j- fsCu^d 



«2,u/uC s /2 + V 2 ,sfsC u 
V2 - s ° = f r /o + f r ' V2 ' n = V2 ' S ' 

Ju^s j & 1 J S^U 



(4) 



with similar expressions for isospin partners and antipar- 
ticles. 

The above relations become simpler if the quanti- 
ties i?2,i, fi, and Ci are independent of isospin and are 
also the same for strange and antistrange quarks, i.e., 
u = d = q, u = d = q, and s — s. These condi- 
tions are approximately satisfied in heavy ion collisions 
at RHIC as the ir + /ir~ ratio is almost one around central 
rapidity ^9) . In this isospin symmetric and strange- 
antistrange symmetric limit, the v 2 values for hadrons at 
a given high px are given by: 



V2,N = V2,q, V 2i N 
V2,ir+ = «2,7r° = V2.7V- 
V 2 , q + r s V2,s 



V2,q, V2 1 <t> = V 2 .n 
V2,q + r<jV 2 ,q 



V2,. 



l+Tq 



V2,K+ 



V 2 ,K~ 



V2,A = V2X = 



1 + r s 

2v2,q + r s V 2 ,s 



V2,q + r s V 2 ,s/rq 

1 + r s /r q 



,J 2,A 



J 2.T. 



2 + r s ' 
2v 2 ,q + r s v 2 , s /rq 



r s /rq 



V2,q/2 + V2, S r s 

1/2 + r s 



V2 : q/2 + V2, s r s /rq 
1/2 + r s /r q ' 



(5) 



with N denoting a nucleon. In the above, the px- 
dependent variables r q and r s arc defined as 



fq c q 



fq c q 



fsCq 
fqC-s 



(6) 



Tq can be determined from 



From Eq.(|5|), we see that, e.g 
the V2 of high px pion, proton and antiproton, while r s 
can be determined from the V2 of high px proton, kaon, 
and <f> meson. 

Since the K + /K~ ratio is close to one and p/p ra- 
tio is about 0.7 in heavy ion collisions at RHIC [E7] H|; 
and they should be closer to one in heavy ion collisions 
at LHC, we consider the case where the variables v%,i, fi 
and Ci are the same for quarks and antiquarks, i.e., q = q 
and thus Tq = 1. For such a quark- antiquark symmetric 
partonic matter, Eq.mft simplifies to: 



f2,7r = 
V2,K = 

V2;. 



V2.N 

V2,q 



V2,q, V2,4> 
r s V2,s 



V2,n 



V2„ 
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''2, 



2r s v 2 „ 



V2,A = V2,Y, 



2v 2 . q + r s v 2 „ 



l + 2r s 



(7) 



(8) 



Eliminating the variable r s in Eq. rtq), we obtain two 
relations involving the «2 values of four different hadron 
species, and they can be any two of the following three 
relations: 
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(v2,tt-V2,k) ■ {V2,A-V 2 ,<f>) = 2 (V2 :7 T-V2A) ' (V2,K- «2,^) , 
2 {v 2 ^-V 2 ,k) ■ (V2,E-V2,<p) = (v 2 ,^2,s) ' (V2,K^"2,4>) , 
(v 2 ,7T-V2,e) ■ (W 2 ,A-W2>) = 4 (W2,7r-"2,A) ' C^2,E ^2,0 ) ■ (9) 

These relations on the elliptic flow of hadrons of differ- 
ent flavors become even simpler in several limits for the 
value of r s . In the limit of r s — > due to f s /f q — > 0, i.e., 
if there are very few high pt strange quarks relative to 
light quarks, we have V2 : n = v%,k = V2,n — «2,A = ua.s = 
«2,h = v 2 , q , v 2 ^ = v 2 m = v 2 , s - This is simply due to the 
fact that all strange hadrons with light valence quarks 
consist of leading light quarks. In the opposite limit of 
r s — > oo, i.e., if the number of high pt strange quarks 
is much larger than that of light quarks or the capture 
probability of a soft strange quark is much smaller than 
that of a soft light quark, all strange hadrons with light 
valence quarks consist of leading strange quarks. In this 
case, we have v 2 ^ = v 2 ,n = v 2 , q , v 2 ,k = v 2 ^ = v 2 ,a = 
v 2 ,t, = v 2 ,s = V2S1 = v 2 , s - 

Another interesting limit is r s = 1, which would be 
the case if the spectrum of high pt strange quarks is the 
same as that of light quarks, and the capture probability 
of a soft strange quark is the same as that of a soft light 
quark, or even though the above two factors are different 
but they cancel each other. In this limit, Eq.(@) gives 



V2,q+02,s 2^2,0+^2,5 
V2,K= ^ -, V 2 \ = 



v 2 .\ 



V2,q+2V 2 ,. 



(10) 
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These relations, together with Eq.(Q), show that the i>2 
values of hadrons at high pt follow a simple quark flavor 
counting rule when r s = 1. 

Eqs. (0-§) show that the dependence of the elliptic 
flows of high px hadrons on their flavor composition is de- 
termined by the relative magnitude of the elliptic flow of 
high pt strange quarks to that of high pt light quarks. If 
strange quarks have the same elliptic flow as light quarks 
at highpT, i.e., t>2, s = V2. q , then v2,h = V2, q for all SU(3) 
hadrons regardless of the value of r s . This is true even if 
the pt spectrum for strange quarks is different from that 
for light quarks. We note that in the present study we are 
only concerned with the relative magnitude of the elliptic 
flow of different hadrons at high pt, not their absolute 
magnitudes or shape. 

On the other hand, fast moving heavy quarks have 
been shown to suffer less energy loss in a thermalized 
parton plasma than fast moving light quarks ]30|-|32]| . In 
the parton transport model, this would imply that high 
Pt strange and heavier quarks may have smaller scatter- 
ing cross sections than light quarks in a partonic matter. 
It is thus possible that the elliptic flow of strange quarks 
is smaller than that of light quarks, i.e., V2, s < ^2,g- In 
this case, we obtain from Eqs. the following flavor 

ordering of the v 2 values for hadrons at a given high pt : 

{v2,tv=V2,n) > {V2,A=V 2 ,Y,) > V 2 ,K>V2,~ > (^2,0=^2,0 ) ■ (H) 

In Fig. [I], we illustrate the flavor ordering of hadron 
elliptic flows at high pt for the case of v 2 . s < «2,g- The 



spacings between different curves correspond to the case 
of r s = 1 and thus follow the quark counting relations of 
Eqs.(0) and ([To|). We note that the vertical scale for v 2 
is in arbitrary units, and the shape of V2 as a function 
of pt is also arbitrary. The scale po denotes the typical 
transverse momentum above which the pt spectra of final 
partons changes from soft to hard, and its value should 
probably be a few GeV/c. All curves are shown well 
above po, reflecting the fact that the relations derived in 
the present study only apply to hadron elliptic flows at 
high pt- In general, r s can take any finite positive value, 
but the flavor ordering of hadron elliptic flows remains 
similar to that shown in Fig. [I] as long as V2, s < «2,g- 
However, the spacings between different curves can be 
different, while still being constraint by the two relations 
given in Eq. (||) . Since the t>2 magnitudes of hadrons fol- 
low the mass ordering at low pt and the flavor ordering 
at high pt, the curve for kaon «2(pt), which is above 
those for proton and A at low pt , will cross and become 
lower than the latter two curves as pt increases. A simi- 
lar relation exists between the curve for <f> meson «2(pt) 
and those for A and S. 
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FIG. 1. Schematic plot for the flavor ordering of the elliptic 
flows of hadrons at high pt- Details are given in the text. 

We have not included the effects of resonance decays 
in this study. The relations shown in Eq.(|J) can be ex- 
tended to resonances such as 77, p, ui, K* , and A. These 
resonances at high pt will decay to stable hadrons at dif- 
ferent transverse momenta, thus complicating the rela- 
tions we have thus derived for hadrons which are directly 
formed from the quark coalescence. Since the transverse 
momentum of a decay product is usually small compared 
to that of the parent hadron, and the inverse power law 
spectrum shows a rapid decrease with pt, we expect that 
the resonance contribution to hadron elliptic flow at high 
Pt is small compared to the contribution from directly 
formed hadrons. 

In summary, using a parton coalescence model to de- 
scribe the formation of hadrons from the initial partonic 
matter in ultra-relativistic heavy ion collisions, we have 
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studied the dependence of the elliptic flows of hadrons 
at high pt on their flavor composition. Since the ellip- 
tic flow is generated mostly in the early partonic phase, 
and high px hadrons are mainly formed from the coales- 
cence of a high pt quark or antiquark produced from the 
initial hard processes and low px quarks or antiquarks 
from the soft processes, the magnitudes of the hadron 
elliptic flows at high px are determined by that of high 
Pt quarks. The relations between hadron and parton 
elliptic flows at high px also depend on the final quark 
spectrum at high px (/i(pr)) and the capture probabil- 
ity of a soft quark (cj) by a high px quark to form a 
high px hadron. If strange quarks have a smaller ellip- 
tic flow than the light quarks, then the quark coales- 
cence model leads to the flavor ordering in the elliptic 
flows of the hadrons formed from an isospin symmet- 
ric and quark-antiquark symmetric partonic matter, i.e., 

(v2,iv=V2,n) > (w 2 ,A=W2,s) > V2.K > V 2 ,B > (l>2 ,<t>= v 2 ,fi) ■ 

We have also obtained two relations which are indepen- 
dent of fi(pr) and Cj and involve the elliptic flows of 
four hadron species at high px- In the special case that 
/ 2 (px) and Ci are the same for strange quarks and light 
quarks, values of the elliptic flows for high px hadrons of 
different flavors are found to follow the quark counting 
rule. It will be very interesting to test these predictions 
in current and future heavy ion collisions. Such studies 
will provide valuable information on whether a partonic 
matter is formed in the collisions and the subsequent for- 
mation of hadrons can be described by the quark coales- 
cence model. 
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Foundation under Grant No. A-1358, and the Texas 
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